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ABSTRACT 

Aims. In this paper we show expectations on the radio-X-ray luminosity correlation of radio halos at 120 MHz. According to the 
"turbulent re-acceleration scenario", low frequency observations are expected to detect a new population of radio halos that, due to 
their ultra-steep spectra, are missed by present observations at ~ GHz frequencies. These radio halos should also be less luminous 
than presently observed halos hosted in clusters with the same X-ray luminosity. 

Methods. Making use of Monte Carlo procedures, we show that the presence of these ultra-steep spectrum halos at 120 MHz causes 
a steepening and a broadening of the correlation between the synchrotron power and the cluster X-ray luminosity with respect to that 
observed at 1.4 GHz. 

Results. We investigate the role of future low frequency radio surveys, and find that the upcoming LOFAR surveys will be able to test 
these expectations. 

Key words. Radiation mechanism: non-thermal - galaxies: clusters: general - radio continuum: general - X-rays: general 



1. Introduction 

Radio halos are diffuse synchrotron sources from the 
intra-cluster medium (ICM) extended on mega-parsec scale 
(e.g., Feretti 2005; Ferrari et al. 2008). They provide the most 
important evidence of non-thermal components (relativistic par- 
ticles and magnetic fields) mixed with the hot ICM. Galaxy 
clusters hosting radio halos are always characterized by a non- 
relaxed dynamical status suggestive of recent or ongoing merger 
events (e.g., Buote 2001; Schuecker et al 2001; Govoni et al. 
2004; Venturi et al. 2008; Giacintucci et al. 2009). Furhermore, 
the halo radio power at 1 .4 GHz increases with the cluster X-ray 
luminosity, mass and temperature (e.g., Liang et al. 2000; EnBlin 
& Rottgering 2002; Bacchi et al. 2003; Clarke 2005; Dolag et al. 
2005; Cassano et al. 2006, 2007; Brunetti et al. 2009; Rudnick & 
Lemmerman 2009; Giovannini et al. 2009). These correlations 
and the radio halo-merger connection suggest that gravity pro- 
vides the reservoir of energy to generate the non-thermal com- 
ponents (e.g., Kempner & Sarazin 2001). Cluster mergers drive 
shocks and turbulence in the ICM that may amplify the magnetic 
fields (e.g., Carilli & Taylor 2002; Dolag et al. 2002; Bruggen 
et al. 2005; Subramanian et al. 2006; Ryu et al. 2008) and ac- 
celerate high energy particles (e.g., Fujita et al. 2003; Hoeft & 
Bruggen 2007; Brunetti & Lazarian 2007; Pfrommer et al. 2008; 
Vazza et al. 2009). 

Two main scenarios have been proposed to explain the ori- 
gin of relativistic particles in radio halos, namely i) the turbu- 
lent re-acceleration model, whereby relativistic electrons are re- 
energized in situ due to the interaction with MHD turbulence 
generated in the ICM during cluster mergers (e.g., Brunetti et 



Send offprint requests to: R.Cassano 
* e-mail: rcassano@ira. inaf. it 



al. 2001; Petrosian et al. 2001), and ii) the secondary electron 
models, whereby the relativistic electrons are secondary prod- 
ucts of the collisions between cosmic rays and thermal protons 
in the ICM (e.g., Dennison 1980; Blasi & Colafrancesco 1999, 
Pfrommer & EnBlin 2004). 

Observations provide support to the idea that turbulence may 
play a role in the particle re-acceleration process (e.g., Brunetti 
et al. 2008; Ferrari et al. 2008; Cassano 2009; Giovannini et al. 
2009), in which case, the population of radio halos is predicted 
to be a mixture of sources with different spectral properties, with 
halos having steeper spectra being more common (Cassano et 
al. 2006a, hereafter C06; Cassano et al. 2009, hereafter C09). 
In this respect, since very steep spectrum halos should glow up 
at low radio frequency, upcoming observations with the Low 
Frequency Array (LOFAR) and the Long Wavelength Array 
(LWA) will be crucial. 

In this paper we discuss how such a predicted population 
is expected to affect the properties of the radio-X-ray lumi- 
nosity correlation at low radio frequency, and investigate the 
potential of LOFAR surveys. A ACDM cosmology (H a = 
70kms-' Mpc' 1 , Q. m = 0.3, Q A = 0.7) is adopted. 



2. The population of ultra-steep spectrum radio 
halos 

The formation and evolution of radio halos according to the tur- 
bulent re-acceleration scenario have been investigated by means 
of Monte Carlo based procedures (Cassano & Brunetti 2005; 
C06; C09). These procedures allow us to account for the main 
ingredients in the model, i.e., the rate of cluster-cluster mergers 
in the Universe, their mass ratio, and the fraction of energy dis- 
sipated during mergers that is channelled into MHD turbulence 
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and acceleration of relativistic particles. We simulate the forma- 
tion history of ~ 1000 galaxy clusters with present day masses 
in the range ~ [0.2 - 6] X 10 15 M@. 

Turbulence acceleration is a rather inefficient process in the 
ICM and electrons can be accelerated up to energy of sev- 
eral GeV, since at higher energy the radiation losses dominate 
(e.g., Brunetti & Lazarian 2007). This implies the presence of 
a gradual spectral steepening at high frequencies in the syn- 
chrotron spectrum of radio halos. The presence of this steepen- 
ing makes it difficult to detect radio halos at frequencies larger 
than the frequency v s at which the steepening becomes severe. 

Following C06 and C09 we use homogeneous models that 
assume: i) an average value of the magnetic field strength in 
the radio halo volume that scales with the cluster mass as B — 
B<m>(M v / < M >) 6 Q ii) that a fraction, r] u of the PdV work, 
done by subclusters crossing the main clusters during mergers 
goes into magneto-acoustic turbulence. The frequency v s , de- 
fined as the frequency at which the spectral slope of the halos 
becomes a > 1.9 (P{v) oc v~ a ), depends on the acceleration ef- 
ficiency in the ICM, on the magnetic field in the ICM, B, and 
on the energy density of the cosmic microwave background radi- 
ation (CMB) as: v s oc (Bx 2 )l(B 2 + B 2 cmb f. The frequency v s is a 
more practical re-definition of the synchrotron break frequency, 
Vb, and, in the case of homogeneous models is v, ~ 7 (C09). 
According to C09, in the case of a single major merger between 
a cluster of mass M v and a subcluster with mass AM, v s is 
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Fig.l. RHLF at 120 MHz (black lines) and in the redshift in- 
terval z = 0-0.1. The differential contributions from halos 
with 120 < v iS < 240 MHz (magenta line, star symbols), 
240 < v s < 600 MHz (green line), 600 < v s < 1400 MHz 
(blue line) and v s > 1400 MHz (red line) are also shown. The 
dashed lines correspond to the minimum detectable halo power 
for £ • F = 0.25, 0.6 and 1 mJy/beam (from left to right, see text). 
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where B cm i, = 3.2(1 + z) 2 //G is the equivalent magnetic field 
strength of the CMB, and R v is the cluster virial radius. It is ex- 
pected that mergers may generate halos with larger v s in more 
massive clusters, and that halos in clusters with the same mass 
M v (and magnetic field) and redshift could have different v s 
depending on the properties of the merger event responsible 
for their generation. Halos with v s > 1 .4 GHz must be gener- 
ated in connection with the most energetic merger-events in the 
Universe since only these mergers may allow for the efficient ac- 
celeration that is necessary to have relativistic electrons emitting 
at these frequencies. Present surveys carried out at v„ ~ 1 GHz 
detect radio halos only in the most massive and merging clus- 
ters (e.g., Buote 2001, Venturi et al. 2008). On the other hand, 
radio halos with smaller values of v s must be more common, 
since they can be generated in connection with less energetic 
phenomena, e.g., major mergers between less massive systems 
or minor mergers in massive systems, that are more frequent in 
the Universe. This has been addressed quantitatively by means 
of Monte Carlo calculations that allow us to derive the fraction 
of clusters with radio halos with different v s as a function of the 
cluster mass and redshift. The expected population of radio halos 
is indeed constituted by a mixture of halos with different spectra, 
with steep spectrum halos being more common in the Universe; 
~ 200 radio halos with 120 < v s < 600 MHz are expected in 
future LOFAR surveys at 120 MHz (C09). 

According to this model the monochromatic luminosity of 
radio halos at a given frequency, vq < v s , increases with increas- 
ing v s . For a fixed cluster mass (or X-ray luminosity), the rela- 



1 B <M> is the value of the magnetic field averaged in a region of ra- 
dius = 500 tcl kpc in a cluster with viral mass < M >= 1.6 x 10 15 M Q . 



tion between the radio luminosity at vo of halos with v s = v s j 
and v Sj 2, P Vtl (vo,L x ) and P Vs2 (vo,L x ) respectively, is (C09) 



P Vsl (v ,L x ) = P Y Jv ,L x )(^X 

\Vs,2> 



(2) 



where a a 1.3 (e.g., Ferrari et al. 2008) is the radio spectral 
index. A correlation between the radio power Pr of halos and the 
mass (and X-ray luminosity) of the hosting clusters is expected 
(C06). In the simplest case that halos are generated by a single 
major merger this is 
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where the parameter F is defined by T oc M r (r 2/3 in the 
virial scaling). Eq. [3] implies that more massive clusters host 
more luminous radio halos. By considering halos with v, > 1.4 
GHz, C06 showed that the slope of this scaling is consistent 
with that of the observed P(1A) - M v correlation, provided that 
the model parameters (B<m>, b, tj,) lie within a fairly constrained 
range of values (see Fig. 7 in C06). We refer the reader to Sects. 
3.3 and 4.1 of C06 for a more detailed discussion on model pa- 
rameters and on their constraints. 

Following C09, we adopt a reference set of parameter^ 
B<m> = 1.9/iG, b = 1.5, rj t = 0.2, that falls in that range and 
sets a M = 3.3, with P(1.4) oc M" M . This implies P(1.4) oc L 2 / 5 
assuming the L x - M v correlation for galaxy clusters as derived 
in C06. Since the bulk of radio halos in our calculations is found 
to be associated with clusters of mass ~ [1 - 2] x 10 15 M Q the 
adopted values of B and b imply typical average magnetic fields 
~ 1 — 3 fiG. These values of B are similar to those derived from 



2 We note that for this particular configuration of parameters even in 
the case of re-acceleration phase the energy of magnetic field is always 
dominant with respect to that of relativistic electrons. 




Fig. 2. Expected distribution of radio halos in the P(120) - Lx diagram (colored open dots, the color code is the same as Fig.Q]) 
together with the position of radio halos observed at 1.4 GHz (black filled dots). The correlation P(1.4) - Lx extrapolated at 120 
MHz (dashed lines) and the best fit of the P(120) - Lx correlation (solid lines) are also shown. From left to right ^ F — 1, 0.6 and 
0.25 mJy/beam. 



rotation measurements (e.g., Govoni & Feretti 2004; Bonafede 
et al. 2010) and equipartition assumption (EnBlin et al. 1998; 
Govoni et al. 2001). 

The observed P( 1 .4) - Lx correlation shows an intrinsic scat- 
ter across the radio luminosity 5P/P - +2 (e.g., Brunetti et al. 
2009). In principle, in our model a scatter in the P(IA) - Lx 
correlation is expected due to the different monochromatic ra- 
dio luminosity of halos with different v s (Eq. [2]). Our calcula- 
tions show that the fraction of clusters hosting radio halos with 
v. s > 3500 MHz is about a few percent, thus we would expect 
SP/P = +1/2 x (3500/1400) 13 ~ 1.7 for halos with v s > 1.4 
GHz, which is in line with the observed scatter. However, we 
stress that there are other possible sources of scatter which are 
difficult to take into account in homogeneous models. These are 
due to e.g., differences in the cosmic rays and magnetic field 
content in clusters with the same mass. 

Once we anchor the luminosity of halos with v s > 1 .4 GHz, 
Pi.4(1.4, L x ), to the observed P(IA) - L x correlation, Monte 
Carlo calculations carried out by considering an observing fre- 
quency vo and Eq.|2]allow us to derive the expected radio halo lu- 
minosity functions (RHLF; see e.g., C06 & C09 for details). As 
an example, Fig. Q] shows the total RHLF obtained from Monte 
Carlo calculations at vo = 120 MHz (black line) and z = — 0.1, 
together with the differential contributions to the RHLF from 
halos with different v s (see figure caption). As expected, radio 
halos with smaller v s mainly contribute to the low-power end of 
the total RHLF, and the peaks of the RHLF of different popula- 
tions move towards low radio powers with decreasing v s . This 
implies that, depending on their sensitivity, surveys at low radio 
frequency will unveil new populations of halos. 

3. Monte Carlo distributions of radio halos in the 
P(120) -L x plane 

The aim of this section is to investigate how the presence of the 
new population of ultra-steep spectrum halos, predicted in deep 
low frequency radio surveys, may affect the radio - X-ray lu- 
minosity correlation of halos at low radio frequency. LOFAR 
will carry out surveys between 15 and 210 MHz in the Northern 
hemisphere with unprecedented sensitivity and spatial resolu- 
tion. Since LOFAR is expected to carry out the deepest large 



area radio surveys at v a = 120 MHz (e.g., Rottgering et al. 2006), 
in this paper we focus on the P(120) - Lx correlation. 

The most crucial point in this respect is the estimate of the 
minimum diffuse flux of a radio halo detectable by these surveys 
as a function of redshift. It is well known that the brightness 
profiles of radio halos smoothly decrease with distance from the 
cluster center (e.g., Govoni et al. 2001; Murgia et al. 2009) im- 
plying that the outermost region of the halos will be difficult to 
detect in radio surveys. Brunetti et al. (2007) found that the typ- 
ical profiles of radio halos are such that about 58% of their flux 
is contained within half radius (Rh)- Following C09, in order to 
have a good sensitivity to diffuse emission, we assume a beam 
of 25 x 25 arcsec and estimate the minimum flux of a Mpc-sized 
halo detectable in a LOFAR survey, fn, by requiring that the 
mean brightness within half radius of a halo, B<o.sr h , is £, times 
the rms (F) of the survey, i.e., 



><0.5 R H 



0.58/ tf 

7T/4 2 H 



f H (z) * 2 x 10" 



{F 



mJy/b 



mJy . 



(4) 



where 8h(z) is the angular size of radio halos, in arcseconds, at 
a given redshift. In the case of £ =s 1 this approach would guar- 
anty the detection, at several <x, of the central brightest region of 
halos, thus leading to the identification of candidate radio halos 
in the survey. 

This simple approach has been tested in Cassano et al. (2008) 
by injecting "fake" radio halos in the (u,v) plane of NVSS and 
GMRT observations. It has been shown that radio halos become 
visible in the images as soon as their flux approaches that ob- 
tained by Eq. |4]with <f = 1 and £ = 2 for the NVSS and GMRT 
observations, respectively (see also Brunetti et al. 2007; Venturi 
et al. 2008). 

The LOFAR Large Area Survey is expected to reach an 
rms ~ 0.1 mJy/beam at 120 MHz with a beam ~ 5x5 arc- 
sec. These resolution and sensitivity are obtained by assuming 
a "uniform weighting", thus we may assume that the same sen- 
sitivity can be reached with a larger beam, obtained for exam- 
ple by e.g., "natural weighting". However, the survey sensitiv- 
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Fig. 3. Radio - X-ray luminosity correlation of giant radio halos at 1.4 GHz. Observed halos at 1.4 GHz (black points) and "sim- 
ulated halos" with v s > 1 .4 GHz (red points) are plotted together with their best-fit relations (solid black line and dashed red line, 
respectively). Left Panel: the observed halos are compared with those expected by the model by considering X-ray and radio sensi- 
tivities, and sky coverage of GMRT+NVSS samples (see text). Right Panel: present observed halos are compared with those, with 
v s > 1 .4 GHz, that are expected to be detected by LOFAR surveys assuming £ F ^ 0.25 mJy/beam. 



ity may be limited by the rms confusion level. This is given by 
(e.g., Condon 1987; Kronberg et al. 2007): 



Tconf - 0.13 X 



#1 X6 2 



25 x 25 arcsec 2 



/ v \-°- 7 

mJy (5) 

\ 120 MHz/ y V ; 



where 6*1,2 is the beam size in arcsec, and v is the frequency in 
MHz. 

Of course all the issues discussed above will be clarified dur- 
ing the commissioning phase of LOFAR. Thus we decided to 
present calculations in several cases, specifically £ • F — 0.25, 
0.6 and 1 mJy/beam to cover a range of possible LOFAR sen- 
sitivities 0. Vertical dashed lines in Fig. Q] show the minimum 
power of a halo at z ~ 0.05 detectable by LOFAR surveys as- 
suming g-F = 0.25, 0.6 and 1 mJy/beam. The important point is 
that with increasing survey sensitivity new populations of radio 
halos are expected to be unveiled, with the detectable number of 
ultra steep spectrum halos increasing in deeper surveys. 

LOFAR observations will allow to study the distribution of 
radio halos in the radio-X-ray luminosity diagram at low ra- 
dio frequencies, so far an unexplored issue. The vast majority 
of ultra-steep spectrum halos visible at low frequencies are ex- 
pected to be associated with galaxy clusters of intermediate X- 
ray luminosity, Lx ~ 3 - 6 ■ 10 44 erg/s, and should be less lu- 
minous than radio halos that are presently observed at GHz fre- 
quencies. This should affect the radio-X-ray luminosity corre- 
lation of halos at low frequencies, that is expected to be steeper 
and with larger scatter than that at 1.4 GHz. 

To address this issue quantitatively we assume vo = 120 
MHz, and following C06 and C09 we use Monte Carlo proce- 
dures based on the extended Press & Schechter (1974; Lacey & 
Cole 1993) formalism to obtain i) the population of galaxy clus- 
ters, with their mass (and X-ray luminosity), in the redshift in- 
terval z = 0-0.5, and ii) the population of radio halos, with their 



3 Note that our choice of £ ■ F is thought to mimic different possible 
configurations, e.g., {; ■ F = 0.6 mJy/beam can be £ = 3 (3<x detection of 
the average halo brightness in half R H ) and F=0.2 mJy/beam, or £ = 1 
and F=0.6 mJy/beam 



v s , associated with these clusters. We use homogeneous models 
and the set of model parameters given in the previous section. 
From these simulations we extract the population of radio halos 
that can be detected by observations at vo = 120 MHz according 
to their radio luminosity and / m ,„(z). 

In particular, the luminosity at 120 MHz of radio halos with 
v s > 1.4 GHz, in clusters with X-ray luminosity Lx, is ob- 
tained from the P(1.4) - Lx correlation, assuming a spectral in- 
dex a = 1.3 and allowing for a random scatter 6P/P = +2 (see 
discussion in Sect.2). The luminosity at 120 MHz of radio halos 
with a given v s is obtained according to Eq. [2] In particular, we 
calculated halo statistics by assuming the following frequency 
ranges: v s = 120 - 240 MHz, 240 - 600 MHz, 600 - 1400 MHz. 
Eq. [2] also implies that halos with v\ < v s < V2 should have 
radio luminosities at 120 MHz which may scatter by a factor 
(vilv\) a that implies 6P/P - +1.3 - 1.7 for the frequency bins 
we are considering. Finally, we assume the LOFAR sky cover- 
age (the Northern hemisphere, 5 > 0, and high Galactic latitudes, 
\b\ > 20) and f min (z) from Eq. g] 

The resulting theoretical distribution of radio halos in the 
P(120) - Lx diagram is shown in Fig. [2 assuming f • F = 1, 0.6 
and 0.25 mJy/beam (colored open dots; from left to right). 
Different colored dots indicate halos with different values of v s 
(the same color code used in Fig.Q]). Halos with different v s fill 
different regions, with radio halos with smaller v s typically lo- 
cated in regions of lower radio luminosities. The number of ha- 
los with smaller v s increases with increasing survey sensitivity. 
In the case of high sensitivity surveys these halos dominate the 
population and their presence affects the overall shape of the cor- 
relation. The correlation at 120 MHz is predicted more scattered 
and steeper than that observed at 1 .4 GHz. A quantitative esti- 
mates of the steepening can be obtained by repeating many times 
the Monte Carlo procedure described above and by fitting the ob- 
tained halo distributions in the P(120)-Lx diagram. Fig.|4]shows 
a histogram of the slopes of the correlation obtained after 100 
Monte Carlo runs in the case tj F ~ 0.25 mJy/beam. The mean 
value of the slope is a corr =t 2.68, while we find a corr ~ 2.45 and 
2.46 in the cases ijF ~ 1 and 0.6 mJy/beam, respectively (with 
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Fig. 4. Spectral slopes of the P(120) - Lx correlation obtained 
after 100 Monte Carlo extractions of the radio halo distribution 
in the P(120) - Lx diagram, assuming £ F — 0.25 mJy/beam. 
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Fig. 5. Fraction of clusters with radio halos with 120 < v s < 240 
MHz (magenta lines) and 600 < v s < 1400 MHz (blue lines) 
as a function of the cluster X-ray luminosity. The calculations 
are showed for the redshift range z = 0.2 - 0.3 and for B<m> = 
1.9 pG, b = 1.5, T], = 0.2 (solid lines) and B <M> = 0.2 pG, 
b = 0.6, T], = 0.38 (dashed lines). 

68% of values typically within Aa ~ 0.07). The values of a corr 
are significantly larger than that at 1 .4 GHz. 

For completeness, in Fig. [3] left panel, we show the theoret- 
ical distribution of giant radio halos with v, > 1400 MHz (red 
points) in the P(1400) -Lx plane together with the observed cor- 
relation of halos at 1400 MHz (black points, taken from Brunetti 
et al. 2009; see Tab. 1 and references therein). In this case, to 
compare model expectations and present observations, we fol- 
low C09 (Sect.3) and derive the theoretical distribution of ra- 
dio halos by considering the combination of the NVSS-XBACs 
(Giovannini et al. 1999) (radio-X-ray) selection criteria and sky 
coverage (at z = 0.044 - 0.2), and the X-ray luminosity range 
and sky coverage of the GMRT radio halo survey (Venturi et al. 
2007, 2008) (at z = 0.2 - 0.32fl The observed and theoreti- 
cal distributions show a good agreement, although data points 

4 The bulk of halos is found in these surveys. 



present a slightly larger scatter than expectations which can be 
easily interpreted as due to variations of magnetic field in clus- 
ters with the same X-ray luminosity (see Sect. 2). In Fig. [3] right 
panel, the same observed distribution of radio halos at 1 .4 GHz 
(black points) is compared with that of "simulated" halos with 
v s > 1400 MHz (red points) detectable by a LOFAR survey at 
120 MHz with £ F=0.25 mJy/beam. As expected, radio halos 
with v s > 1 .4 GHz follow a trend consistent with presently ob- 
served halos, while their larger number simply reflects the large 
sensitivity of LOFAR surveys with respect to present surveys. 

3.1. Dependence on model parameters 

The steepening of the correlation is independent of the adopted 
values of model parameters, at least by considering sets of pa- 
rameters in the region (B < m>, b, rj t ) that reproduce both the ob- 
served slope of the P(1.4) - M v correlation (a M = 2.9 ± 0.4) 
and the observed fraction of galaxy clusters with radio halos. 
C06 and C09 already discussed the dependence of expectations 
on model parameters. They showed that the expected number of 
radio halos decreases only by a factor of ~ 2 - 2.5, from super- 
linear (b > 1) to sub-linear (b < 1) magnetic scaling (see also 
Fig. 4 in Cassano et al. 2006b). 

For a fixed value of b larger values of B<m> produce P( 1 .4) oc 
M" M correlations with slightly flatter slopes (see Tab.3 in C06). 
For example, in the case b = 1.5 the allowed values of B range 
from B<m> — 1 -9 pG to ^ 2.8 pG, and correspondingly it is «m = 
3.3 and 2.5, respectively, still consistent with the observed one. 

The steepening of the correlation is due to the glow up of 
new radio halos at low frequency, thus another point is whether 
the fraction of halos with smaller v s (ultra steep spectrum halos) 
changes from super-linear to sub-linear cases. To investigate this 
effect, in Fig. [5] we report the percentage of radio halos with 
120 < v., < 240 MHz (magenta lines) and 600 < v s < 1400 
MHz (blue lines) as a function of the cluster X-ray luminos- 
ity. In Fig. [5] we assume two configurations of parameters : the 
one used in the present paper (solid lines), and a sub-linear one 
(B <M> = 0.2 pG, b = 0.6, T], = 0.38, dashed lines). In both 
cases, the vast majority of radio halos hosted in clusters with 
L x < 3 x 10 44 erg/s has v s < 240 MHz, while halos with v s > 600 
MHz become dominant in more luminous clusters. On the other 
hand, we find that the fraction of halos with 120 < v s < 240 
MHz is larger in the sub-linear case. This is because it is more 
difficult to generate radio halos with larger v s in the case of lower 
magnetic fields (provided that radiative losses are dominated by 
the Inverse Compton losses due to the CMB photons). We may 
conclude that in sub-linear cases we expect the following main 
effects: i) the f(120) - Lx diagram should be less populated 
than that in super-linear cases (less halos are expected); ii) the 
P( 1 .4) - Lx correlation is expected to be flatter than that in super- 
linear cases (see Tab. 3 in C06); Hi) the P(120) - Lx diagram 
should be even more dominated by ultra-steep spectrum halos, 
making the steepening of the correlation at lower frequency even 
stronger. 

4. The effect of an evolving magnetic field 

In this paper we have applied a statistical model based on the 
turbulence acceleration scenario (discussed and developed in 
Cassano & Brunetti 2005, C06 and C09) to derive the expected 
distribution of radio halos in the P(120) - Lx diagram. The cos- 
mological evolution of the magnetic field in this model is ac- 
counted for by scaling the field with the cluster mass, as sug- 
gested by cosmological MHD simulations (e.g., Dolag et al. 
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Fig. 6. Left Panel time-evolution of the ratio between v s computed by assuming an evolving (as in Fig. 2 of SSH06) and a constant 
magnetic field; Central Panel time-evolution of the synchrotron spectra computed by assuming an evolving (dashed lines) and a 
constant (solid lines) magnetic field. Spectra are reported for t = 0.9, 1.1, 1.5, 1.8 and 2 Gyr (see color code in the figure panel). 
Right Panel synchrotron spectral indices between 1.4 GHz and 330 MHz of the synchrotron spectra reported above: for a constant 
magnetic field (solid line) and for an evolving magnetic field (dashed line). 



2002; see also Sect.2). On the other hand, our calculations do 
not follow self-consistently turbulence and the amplification of 
magnetic fields due to this turbulence. The main reason for that 
is that in the turbulence acceleration scenario radio halos are 
generated and disappear due to the acceleration and cooling 
of the emitting particles, and such processes are much faster 
(~ 0.1- 0.2 Gyr) than the slow decay of the magnetic field in the 
ICM, ~ few Gyr (e.g., Brunetti et al. 2009). Moreover, although 
particle-acceleration is demonstrated to be connected with clus- 
ter mergers {e.g., Buote 2001; Govoni et al. 2004; Venturi et al. 
2008), several mechanisms/sources of magnetic field, other than 
mergers-induced amplification, may significantly contribute to 
the magnetic field in the ICM. As a matter of fact stronger values 
of magnetic fields are measured in cooling core clusters, that are 
not merging systems {e.g., Carilli & Taylor 2002; Govoni 2006). 
Although it is clear that a self-consistent treatment of turbulence, 
particle acceleration and magnetic field evolution is mandatory 
and deserves future theoretical efforts, in this Section we show 
that the results presented in this paper are not substantially af- 
fected by the evolution of the magnetic field when clusters be- 
come more relaxed after a merging phase. 

In order to investigate the effect of an evolving magnetic field 
on our results, we consider the simulations of the generation and 
decay of dynamo-active turbulence developed by Subramanian, 
Shukurov & Haugen (2006, hereafter SSH06). They studied the 
decay phase of an induced turbulent flow and magnetic field, that 
follows a saturation phase after the driving is switched off. They 
found that after the exponential growth of the magnetic field a 
saturation phase follows, then the turbulent energy and the mag- 
netic energy decay in such a way that after an eddy turnover 
timescale the turbulent energy density is « 1/2 of its satura- 
tion value, while the magnetic field strength is still at « 90% of 
the saturation value (see Fig. 2 in SSH06). By considering ran- 
dom motion with a typical initial speed vo = 500 km/s and scale 
Iq 300 kpc, which are appropriate for cluster-merger driven 
turbulence, the eddy turnover timescale is to — 0.6 Gyr. Since 
we are interested in the evolution of radio halos on timescales 
of ~ Gyr, and the radiative lifetime of the emitting particle is 
of ~ 0.1 - 0.3 Gyr, we can estimate the variation of the fre- 
quency v 5 as a function of time by adopting the classical formula 



v s oc {Bx 2 )l{B 2 + B 2 cmb ) 2 (see also Sect. 2). The ratio between v s 
in the case of time-dependent magnetic field and in the case of 
constant magnetic field is reported in Fig. [6] (left panel), where 
to = 1.1 Gyr (since here we have considered a saturation phase 
of At = 0.5 Gyr before the driving is switched off). 
Substantial differences are found only for t > 2 Gyr when, how- 
ever, the halo is expected to have already disappeared (being 
v s oc ef, with e, being the turbulence energy density). The ratio 
between the steepening frequencies (Fig. [6] left panel) at t ~ to 
gives an estimate of the error we can make by neglecting the 
time-evolution of the magnetic field, that is < 30%. 

In the central panel of Fig. [6] we also report the time- 
evolution of the synchrotron spectra at different time (? = 0.9, 
1.1, 1.5, 1.8 and 2 Gyr, see figure caption) by assuming a con- 
stant magnetic field (as in the adopted model, solid lines) and an 
evolving magnetic field (as in SSH06, dashed lines). For f < 1.1 
Gyr the difference in the monochromatic radio luminosities is 
less then 10%. In addition, the right panel of Fig. [6] shows the 
evolution of the spectral index between 1.4 GHz and 330 MHz 
of these synchrotron spectra. Up to 1.2-1.3 Gyr, the difference in 
the spectral indeces is very small, less then 10%. 

5. Conclusions 

The observed correlations between the halo radio power (at 1 .4 
GHz) and the cluster X-ray luminosity, mass and temperature, 
and the observed connection between radio halos and cluster 
mergers, suggest a link between the gravitational process of clus- 
ter formation and the generation of radio halos. Radio halos are 
likely generated during cluster-cluster mergers where a fraction 
of the gravitational energy dissipated is channelled into the ac- 
celeration of relativistic particles. A crucial expectation of the 
turbulent re-acceleration scenario, put forward to explain ra- 
dio halos, is that the synchrotron spectrum of halos is charac- 
terized by a cut-off at frequency v > v s with v s determined by 
the efficiency of the acceleration process. The presence of this 
cut-off causes a bias, so that present radio observations at ~ 
GHz frequencies are expected to detect only the most efficient 
radio phenomena in clusters, leaving unexplored a large popu- 
lation of radio halos characterized by spectral cut-off at lower 
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frequencies (e.g., C06, Brunetti et al. 2008; C09). Future low 
frequency radiotelescopes as LOFAR and LWA are expected to 
unveil the populations of ultra steep spectrum radio halos, with 
v 5 < 1 GHz, in clusters, providing to test the idea of turbulent 
re-acceleration. 

One may wonder whether halos with v, < 1 GHz could 
be detected by present radiotelescopes at 1 .4 GHz. To address 
this point in Fig. [7] we report the flux distribution at 1.4 GHz 
of halos with 600 < v s < 1400 MHz (that in our model have 
a « 1.7 between 120 MHz and 1400 MHz) that are expected 
to be detected by LOFAR at 120 MHz assuming £F ^ 0.25 
mJy/beam. Calculations are derived by assuming: i) at z < 0.3 
the X-ray flux limit and sky coverage of the extended ROSAT 
Brightest Cluster Sample (eBCS, Ebeling et al. 1998, 2000) and 
of the ROSAT-ESO Flux Limited X-ray Galaxy Cluster Survey 
(REFLEX, Boringher et al. 2004), and ii) at z = 0.3 - 0.6 the 
X-ray flux limit and sky coverage of the Massive Cluster Survey 
(MACS, Ebeling et al. 2001). Calculations show that potentially 
very deep pointed observations at 1 .4 GHz of all these clusters 
may lead to the detection of a few of these ultra steep spectrum 
halos (those with 600 < v s < 1400 MHz); the ultra steep spec- 
trum halo detected in the cluster Abell 521 (Brunetti et al. 2008; 
Dallacasa et al. 2009) belong to this class of halos, although it is 
among the flatter spectrum objects in this class, with v s « 1200 
MHz. 

In this paper, we discuss the consequence of this new popu- 
lation of radio halos on the slope of the radio-X-ray luminosity 
correlation at low frequency. According to homogeneous mod- 
els, ultra-steep spectrum halos are expected to be less luminous 
than halos with larger v s associated with clusters of the same 
mass. Also, radio halos with smaller v s should be statistically 
generated in clusters with smaller mass (and Lx). The combi- 
nation of these two expectations implies that the radio- X-ray 
luminosity correlation should be broader and steeper at lower 
frequencies. 

Based on this model, we perform Monte Carlo simulations 
of the distribution of radio halos in the f(120) - Lx plane. We 
find that halos distribute in the P(120) - Lx plane according to 
a correlation which is steeper (Aa « 0.4) and broader than that 
observed at 1 .4 GHz, with ultra-steep spectrum halos broaden- 
ing the scatter in the region of low luminosity. We find that the 
number of ultra-steep spectrum halos increases with increasing 
the survey sensitivity and this further steepens the correlation. 
The forthcoming LOFAR surveys should constrain the expected 
steepening of the correlation and test our expectations. 

Although a self-consistent treatment of turbulence accelera- 
tion and amplification of the magnetic field in clusters is manda- 
tory and deserve future efforts, we show that the main ingredi- 
ents in the adopted scenario are two "fast" processes: particle 
acceleration and particle cooling that follow the decay of turbu- 
lence. Being a "slow" process, we show that the possible decay 
of the field with turbulence is not expected to affect the modeling 
of halo statistics significantly. 
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